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group  of  atoms  as  a function  of  its  composition,  geometry,  and  local  chemical 
environment.  In  transition  metals  and  transition-metal  coordination 
complexes,  where  local  magnetic  spin  polarization  of  electrons  is  important, 
the  concept  of  orbital  electronegativity  can  be  further  generalized  to  the 
individual  spin-orbitals.  By  viewing  a transition-metal  surface,  cluster,  or 
coordination  complex  as  providing  orbital  or  spin-orbital  pathways  for 
electrons  to  effectively  flow  between  reactants,  where  such  flow  directly 
between  reactants  in  the  gas  phase  is  forbidden  by  orbi tal -symmetry  restric- 
tions or  unfavorable  electronegativity  differences,  orbital  or  spin-orbital 
electronegativity  can  be  used  in  conjunction  with  -SCf=fer  calculations  for 
representative  clusters  and  complexes  as  an  approximate  index  of  hetero- 
geneous or  homogeneous  reactivity.  Recent  applications  of  this  concept  to  a 
number  of  problems  associated  with  reactivity  at  transition-metal  interfaces 
are  reviewed,  including: 

nr  ^e  dissociation  and  reactivity  of  hydrogen  at  low-coordination 
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transition-metal 


sitesj 


The  interaction  of  atomic  hydrogen  with  transition-metal  interfaces^ 
Jhe  surface  reactivity  of  iron 


I 


X-a.\ph<3 


security  Classification  of  this  PAocrw***  Dm  Emwij 


> 

I.  INTRODUCTION 

One  of  the  major  goals  of  surface  chemistry  and  physics  is  to  elucidate 
the  fundamental  nature  of  chemical  reactivity  at  surfaces  and  interfaces. 

In  particular,  transition  metals  and  their  alloys  have  been  a focus  of 
attention  in  regard  to  the  generally  high  surface  reactivity  of  these 
materials,  in  both  their  crystalline  form  and  in  the  form  of  small  particles 
and  clusters  which  constitute  the  active  centers  of  supported  heterogeneous 
catalysts.  Molecular  transition-metal  coordination  complexes  are  often 
centers  of  homogeneous  catalysis. 

In  characterizing  surface  reactivity  and  heterogeneous  catalytic 
activity,  it  is  customary  for  one  to  distinguish  between  those  reactions 
which  are  structure-sensitive  or  "demanding"  and  those  reactions  which  are 
structure-insensitive  or  "facile"  [1]  Most  heterogeneous  reactions  are,  in 
fact,  facile.  The  relatively  few  which  are  demanding  usually  vary  in 
activity  by  no  more  than  one  order  of  magnitude  over  a range  of  surface 
structures,  for  fixed  average  surface  composition.  On  the  other  hand, 
heterogeneous  reactivity,  facile  or  demanding,  may  vary  by  several  orders 
of  magnitude  with  changes  in  surface  composition.  For  example,  among  the 
Group-VIII  transition  metals,  osmium,  iridium,  and  platinum,  the  catalytic 
activity  for  ethane  hydrogenolysis , a demanding  reaction,  varies  by  seven 
orders  of  magnitude  [2].  As  another  striking  example,  the  alloying  of  only 
five  atomic  percent  of  copper  with  nickel  reduces  the  catalytic  activity 
of  the  latter  metal  for  ethane  hydrogenolysis  by  three  orders  of  magnitude, 
attributed  largely  to  surface  segregation  of  copper  [2]. 

The  dominance  of  surface  composition  over  surface  geometry  in  determining 
heterogeneous  reactivity  suggests  that  one  might  look  for  electronic  indices 
of  surface  reactivity,  dependent  on  surface  electronic  structure,  which 
would  be  an  approximate  gauge  of  the  relative  activities  and  selectivities 
of  surfaces  of  different  composition,  with  respect  to  specific  reactants. 
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The  establishment  of  such  indices  could  ultimately  serve  as  a guide  in  the 
systematic  optimization  of  surface  activity  and  selectivity  through  alloying 
or  chemical  modification.  Since  surface  electronic  structure  can  vary  with 
morphology,  e.g.,  surface  "roughness"  or  particle  dispersion,  the  structure 
sensitivity  of  certain  reactions,  although  more  subtle  than  composition 
sensitivity,  should  also  be  within  the  scope  of  such  electronic  indices. 

During  the  past  two  years,  we  have  been  investigating  the  electronic 
structures  of  clusters  and  coordination  complexes  which  are  theoretical 
models  for  active  sites,  chemisorption,  and  reaction  intermediates  at  transition- 
metal  interfaces,  utilizing  the  self-consistent-field  X-alpha  scattered-wave 
(SCF-Xa-SW)  density-functional  approach  to  molecular-orbital  theory  [3],  The 
assumptions  underlying  this  approach  are  that  chemisorption  and  heterogeneous 
reactivity  on  the  active  sites  of  transition-metal  clusters  and  surfaces  are 
governed  by  essentially  the  same  types  of  electronic  factors  which  determine 
the  metal-ligand  bonding  and  homogeneous  reactivity  of  isolated  transition- 
metal  coordination  complexes.  The  initial  phases  of  this  work  are  described 
in  detail  in  two  recent  articles  [4,5]. 

II.  SPIN-ORBITAL  ELECTRONEGATIVITY  AND  THE  SCF-Xa  METHOD 

A recent  outgrowth  of  these  theoretical  studies  is  the  establishment 
of  "spin-orbital  electronegativity,"  defined  by  the  SCF-Xa  spin-orbital 
energies  for  representative  surface  clusters  and  coordination  complexes, 
as  a reactivity  index  of  the  type  described  in  the  preceding  section.  The 
concept  of  spin-orbital  electronegativity  is  derived  from  the  fact  that  the 
orbital  energy  eigenvalues  in  the  SCF-Xa  theory  are  rigorously  equal  to 
first  derivatives  of  the  total  energy  with  respect  to  orbital  occupation 
number  [3,6],  i.e.. 


eiXa  * 3<EXc>/3"i 


(1) 


These  quantities  should  not  be  identified  with  the  orbital  energies  defined 
in  conventional  Hartree-Fock  theory  [6],  namely,  as  the  differences 

ciHF  = <EHF^ni=1^>  " <EHF^ni=0^>  ^ 

between  single-determinant  total  energies  calculated  when  the  ith  orbital  is 

occupied  and  when  it  is  empty  (fixing  the  remaining  occupied  orbitals). 

The  Xa  orbital  energies  defined  in  Eq.  (1)  correspond  closely  to  the  orbital 

electronegativities 


Xi  = 9E/3ni  (3) 

defined  by  Hinze  et  al . [7]  as  a generalization  of  Mulliken's  [8]  definition 
of  electronegativity 

Xm=1<I+A),  (4) 

where  I is  the  ionization  potential  and  A is  the  electron  affinity  of  a 
chemically  bonded  atom  in  its  valence  state. 

This  generalization,  its  relationship  to  SCF-Xa  theory,  and  ultimately 
its  use  as  an  index  of  reactivity  can  be  understood  better  if  one  recalls 
that  electronegativity,  as  originally  defined  by  Pauling  [9],  is  a measure 
of  the  power  of  a chemically  bonded  atom  to  attract  electrons  to  itself. 

Pauling  believed  that  electronegativity  is  a virtually  constant  atomic 
property,  even  for  different  valence  states  of  the  same  element,  and  established 
a scale  of  electronegativities  for  the  elements  based  on  the  empirical  bond 
energies  of  heteronuclear  diatomic  m"',''^ules.  Despite  the  arbitrariness  of 
this  scale  and  the  uncertainties  in  the  thermochemical  data  on  which  it  is 
based,  a wide  variety  of  chemical  phenomena  have  been  reasonably  explained 
through  the  use  of  Pauling's  electronegativity  scale. 


Pauling's  concept  of  electronegativity  as  a fixed  atomic  characteristic 
is  somewhat  more  restrictive  than  Mulliken's  definition  of  electronegativity 
in  terms  of  the  valence-state  ionization  potential  I and  electron  affinity  A, 
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since  one  does  not  expect  that  I and  A or  their  average  should  be  constant 
for  different  valence  or  oxidation  states  of  a chemically  bonded  atom.  Thus 
while  Mulliken’s  electronegativity  scale,  in  its  simplest  form,  can  be 
adjusted  to  agree  reasonably  well  with  Pauling's  scale,  element  by  element, 
Mulliken's  concept  is  more  satisfying  from  a theoretical  point  of  view  and 
allows,  in  principle,  for  the  dependence  of  electronegativity  on  the  chemical 
environment  of  an  atom. 

Since  I and  A are  quantities  related  respectively  to  the  removal  of  an 
electron  from  the  highest  occupied  atomic  orbital  and  the  addition  of  an 
electron  to  the  lowest  unoccupied  orbital,  it  might  be  expected  that 
Mulliken's  concept  of  electronegativity  could  be  further  generalized  to  all 
the  orbitals  of  a chemically  bonded  atom  and  indeed  to  the  molecular  orbitals 
of  an  aggregate  of  atoms.  Thus  one  is  led  to  the  concept  of  orbital 
electronegativity  as  a measure  of  the  power  of  a chemically  bonded  atom  or 
molecular  aggregate  to  attract  an  electron  to  a particular  atomic  or  molec- 
ular orbital.  The  mathematical  definition  of  orbital  electronegativity  as 
the  first  derivative  of  the  total  energy  with  respect  to  occupation  number, 
given  in  Eq.  (3)  in  the  form  suggested  by  Hinze  et  al . [7],  is  consistent 
with  the  above  conceptual  definition.  Implicit  in  Eq.  (3)  are  the  two 
assumptions:  (a)  that  the  occupation  numbers  n^  may  have  both  integral  and 

non-integral  values,  and  (b)  that  once  assumption  (a)  is  made,  then  the 
total  energy  E is  a continuous  and  differentiable  function  of  the  occupation 
numbers . 

The  occupation  numbers  n.  and  statistical  total  energy  E^defined  in 
the  Xa  density-functional  self-consistent-field  theory  fulfill  both  of  the 
above  conditions,  so  that  one  can  uniquely  identify  the  SCF-Xa  electronic 
energy  eigenvalues  EiXa  of  an  atom,  molecule,  or  cluster,  as  given  in 
exhression  (1),  with  the  orbital  electronegativities  defined  in  Eq.  (3). 
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In  the  limit  where  the  total  energy  is  a quadratic  (parabolic)  function  of 
the  occupation  number,  the  orbital  electronegativity  reduces  exactly  to 
Mulliken's  definition  of  electronegativity  given  in  Eq.  (4)  [7].  This 
follows  from  a simple  geometric  theorem  which  relates  the  slope  of  the  chord 
of  a parabola  to  the  slope  of  the  parabola  at  its  midpoint.  The  same  type 
of  argument  applied  to  SCF-Xa  orbitals  leads  to  Slater's  transition-state 
concept  [6],  whereby  one  determines  I or  A for  an  atom,  molecule,  or  cluster 
by  subtracting  or  adding  one-half  a unit  of  valence  orbital  electronic 
charge  and  then  calculating  self-consistently  the  energy  of  the  relaxed 
orbital.  While  these  relaxed  transition-state  energies  can  be  individually 
identified  with  the  corresponding  orbital  ionization  potentials  or 
electron  affinities,  the  unrelaxed  ground-state  SCF-Xa  orbital  energies 

define  a set  of  orbital  electronegativities.  Thus  the  relative  positions 
of  the  SCF-Xa  electronic  energy  levels  for  a system  of  interacting  or  reacting 
atoms,  molecules,  or  clusters  are  a measure  of  the  orbital  electronegativity 
and  chemical -potential  differences  between  the  various  reactants. 

In  those  systems  where  magnetic  spin  polarization  is  important,  one  may 
use  the  spin-unrestricted  version  of  the  SCF-Xa  method  to  calculate  different 
orbitals  for  different  spins,  leading  to  spin-polarized  energy  levels  e^t 
and  elXak  [3,6].  If  these  spin-dependent  orbital  energies  are  identified 
with  orbital  electronegativities,  then  one  is  automatically  led  to  the  concept 
of  spin-orbital  electronegativity  as  a measure  of  the  power  of  an  atom  or 
molecular  aggregate  to  attract  an  electron  to  a particular  atomic  or  molecular 
spin  orbital.  For  example,  the  spin  dependence  of  orbital  electronegativity 
is  central  to  understanding  the  surface  reactivity  of  iron  (see  Section  V). 

To  understand  the  relationship  between  spin-orbital  electronegativity. 


as  represented  by  the  SCF-Xa  spin-orbital  energy  levels  of  representative 
surface  clusters  and  reactants,  and  surface  reactivity,  it  is  helpful  to 
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recdll  the  following  concept  originally  introduced  by  Fukui  [10].  For  a 
concerted  chemical  reaction  to  occur  with  reasonable  activation  energy, 
electrons  must  be  able  to  flow  between  the  reactants  from  occupied  orbitals 
into  unoccupied  orbitals  with  which  they  have  net  positive  overlap,  as  the 
reactants  move  along  the  reaction  coordinate.  Overlap  and  electron  flow 
will  be  ensured  if  the  pertinent  reactant  orbitals  have  the  following 
characteristics:  (a)  the  same  symmetry  (i.e.,  orbital  symmetry  conservation 

as  originally  proposed  by  Woodward  and  Hoffmann  [11]),  and  (b)  equal  or 
nearly  equal  orbital  electronegativities.  In  the  limit  where  electron  flow 
between  reactants  is  simply  from  the  highest  occupied  molecular  orbital 
(HOMO)  to  the  lowest  unoccupied  molecular  orbital  (LUMO),  condition  (b)  is 
equivalent  to  the  requirement  that  the  energy  gap  between  HOMO  and  LUMO  be 
as  small  as  possible  [12]. 

Where  direct  electron  flow  between  reactants  is  forbidden  by  orbital 
symmetry  restrictions  or  unfavorable  orbital  electronegativity  differences 
(implying  a large  activation  energy),  a surface  can  heterogeneously  catalyze 
the  reaction  by  providing  a pathway  for  such  electron  flow,  e.g.,  through 
chemisorption  via  spatially  directed  or  hybridized  d-orbitals  in  the  case 
of  a transition-metal  surface.  Similar  arguments  are  applicable  to  isolated 
transition-metal  coordination  complexes  which  homogeneously  catalyze  electron 
flow  between  reactants  through  bonding  and  exchange  of  ligands. 

III.  THE  DISSOCIATION  AND  REACTIVITY  OF  HYDROGEN 
AT  LOW-COORDINATION  TRANSITION-METAL  SITES 

It  is  well  known,  for  example,  that  certain  "coordinatively  unsaturated" 
transi tion-metal  complexes  in  solution  can  homogeneously  catalyze  chemical 
reactions  [13],  while  it  has  long  been  suspected  that  low-coordination  sites 
on  transition-metal  surfaces  and  supported  transition-metal  clusters  are 
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centers  of  heterogeneous  reactivity  [14,15].  In  this  section,  it  will  be  shown 
that  the  SCF-Xa  electronic  structure  of  such  complexes,  in  conjunction  with  the 
concept  of  orbital  electronegativity,  is  consistent  with  their  reactivity  and 
is  suggestive  of  how  low-coordination  sites  on  transi tion-metal  surfaces  can 
act  as  centers  of  reactivity.  The  dissociation  and  reactivity  of  H2  is 
considered  as  an  illustrative  example. 

As  a working  model,  we  consider  a Group-VIII  transi tion-metal  atom  (M) 
dihedral ly  coordinated  by  ligands  (L),  yielding  the  coordinati vely  unsaturated 
l_2M  complex  illustrated  at  the  top  of  Fig.  1.  This  model  has  the  advantage 
that  it  can  realistically  represent  transition-metal  complexes  of  the  type 
(e.g.,  M = Pt,  Ir,  Rh;  L = Ph^P  = triphenylphosphine)  which  dissociati vely 
bind  and  homogeneously  catalyze  reactions  of  H2  [13],  and  it  can  simulate 
low-coordination  sites  (e.g.,  "corner  atoms")  of  faceted  transition-metal 
clusters  or  stepped  transition-metal  surfaces  which  dissociati vely  chemisorb 
and  heterogeneously  catalyze  reactions  of  H2  [14,15].  In  the  latter  systems, 
the  ligand  (L)  is  also  a metal  atom,  either  of  the  same  species  as  the 
transition  metal  (M),  or  of  a different  species  in  the  case  of  an  alloy  surface 
or  bimetallic  cluster. 

Molecular-orbital  calculations  have  been  carried  out  for  l^M  and  L2MH2 
complexes  by  the  SCF-Xa-SW  method  as  a function  of  metal  species  (M  = Pt,  Ir), 
ligand  species  (L  = phosphine,  Pt),  and  molecular  geometry.  The  resulting 
orbital  energies  for  M = Pt,  L = phosphine,  and  geometry  characteristic  of  the 
platinum-phosphine  complexes  described  in  Ref.  13  are  shown  in  Fig.  1.  Also 
shown,  for  comparison,  are  the  SCF-Xa  orbital  energies  for  the  isolated  metal, 
ligand,  and  hydrogen  molecule  at  the  free-molecule  internuclear  distance 

O O 

0.74  A (H^)  and internuclear  distance  2.8  A (H2*)characteristic  of  the  partially 
dissociated  ("dihydride")  configuration  of  H2  in  the  L2MH2  complex. 

The  SCF-Xa  orbital  energy  eigenvalues  shown  in  Fig.  1 can  be  rigorously 
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identified  with  orbital  electronegativities  which  are  a measure  of  the 
relative  average  electron  donor-acceptor  character  of  the  individual 
orbitals,  as  described  in  the  preceding  section.  Thus  the  fact  that  the 
isolated  ligand  energy  level,  which  corresponds  to  a phosphine  "lone-pair" 
orbital,  nearly  coincides  with  the  d-orbital  energy  level  of  the  isolated 
Pt  atom  (neglecting  relativistic  shifts)  implies  a predominantly  covalent 
L-Pt(5d)  interaction  similar  to  that  expected  for  a direct  Pt(5d)-Pt(5d) 
interaction.  In  this  respect,  the  effect  of  coordinati vely  unsaturated 
phosphine  ligands  on  the  electronic  structure  of  a platinum  atom  is  expected 
to  be  similar  to  that  of  embedding  a Pt  atom  in  a low-coordination  Pt 
environment,  such  as  that  provided  by  a surface  or  cluster. 

The  ligand-metal  interaction  in  the  L2M  complex  leads  to  the  bonding 

orbital  energies  labeled  L-M(d^z)  and  L-M(dz2)  in  Fig.  1,  and  to  the 

antibonding  orbital  energies  labeled  M(dz2)-L*,  M(d^z)-L*,  and  M(s)-L*, 

of  which  M(dyZ)-L*  is  the  highest  occupied  energy  level  in  the  ground  state 

of  the  complex.  A simple  interpretation  of  the  position  of  the  latter  energy 

level  is  that  the  strong  ligand-field  repulsion  of  the  metal  d-orbital  pointed 

on  the  ligand  directions  (the  dyz  orbital  for  the  chosen  coordinate  system) 

raises  the  energy  level  of  this  orbital,  reduces  the  corresponding  orbital 

electronegativity,  and  mixes  in  significant  antibonding  ligand  character. 

The  dz2  orbital  is  also  subject  to  some  antibonding  ligand-field  repulsion, 

whereas  the  do  2>  d , and  dyu  orbitals  remain  essentially  nonbonding, 
x xz  xjr 

When  platinum  atoms  are  substituted  for  the  phosphine  ligands,  the  electronic 
structure  reduces  to  the  manifold  of  bonding,  nonbonding,  and  antibonding 
d-orbital  energy  levels  (the  "d-band")  characteristic  of  a small  platinum 
cluster  [4].  In  this  case,  the  M( dyz ) -L*  (L  = M)  antibonding  orbital  may  be 
interpreted  as  the  analogue  of  a localized  "surface  state"  which  is  split  off 
from  the  top  of  the  d-band. 
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The  most  important  result  of  the  strong  ligand-metal  antibonding 
component  is  to  bring  the  M(d^z)-L*  orbital,  the  highest  occupied  orbital, 
closer  in  energy  and  electronegativity  (as  compared  with  the  isolated  Pt 
atom)  to  the  empty  antibonding  ou  orbital  of  the  molecule.  This  facilitates 
overlap  and  electron  flow  between  the  M(d  )-L*  and  ou  orbitals,  which  are 
symmetry  conserving  [11],  thereby  promoting  dissociation  of  H^.  The  partially 
dissociated  molecule  (H^*),  characterized  by  og  and  ou  orbital  energies  ap- 
proaching the  SCF-Xa  Is  orbital  energy  of  a free  hydrogen  atom  (see  Fig.  1), 
can  bind  or  "chemisorb"  in  a dihydride  configuration  to  the  coordinati vely 
unsaturated  metal  site.  This  is  revealed  by  the  L^MH^  molecular-orbital 
energies  shown  in  Fig.  1 and  the  corresponding  orbital  wavefunction  contour 
maps  shown  in  Fig.  2.  The  2b^  orbital,  for  example,  results  from  overlap 
and  electron  flow  between  the  M(dyZ)-L*  orbital  and  the  ou  orbital.  The 
dihydride  conf iguration  is  further  stabilized  by  the  "butterfly-1  ike"  la^ 
and  2a-j  orbitals  shown  in  Fig.  2,  formed  from  the  overlap  of  the  equatorial 
parts  of  the  L -M( dz2 ) and  M(dz2)-L*  orbitals  with  the  H ( 1 s ) (or  H2*  og) 
orbitals.  Note  that  the  M(dz2)  lobe  pointed  along  the  z-direction  acts  as 
a repulsive  barrier  which  helps  to  keep  the  H atoms  apart.  These  dihydride 
bonding  orbitals  are  offset  somewhat  by  the  4a^  and  3a1  orbitals  resulting 

respectively  from  the  antibonding  interaction  of  the  L-M(d  o)  and  M(d  2 2) 

z x -y 

, orbitals  with  the  H^*  og  orbital,  as  is  evident  in  the  4a^  orbital  contour 

map  shown  in  Fig.  2.  There  is  negligible  contribution  of  the  M(s)  orbital 
component  in  the  binding  of  hydrogen  to  these  platinum  and  iridium  complexes. 
This  is  consistent  with  the  finding,  based  on  SCF-Xa  cluster  calculations  [16] 
and  photoemission  studies  [17],  that  the  metal  d-orbitals  are  almost  exclusively 
responsible  for  the  chemisorption  of  hydrogen  on  second-  and  third-row 
transition  metals  such  as  palladium  and  platinum,  whereas  significant  metal 
s ,d-hybridization  (with  the  s-orbital  component  dominant)  is  involved  in 
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hydrogen  chemisorption  on  first-row  transition  metals  such  as  nickel.  Since 
the  deuterium  molecule  ( ) is  chemically  identical  to  the  hydrogen  molecule 
(H2),  all  the  results  described  above  for  the  dissociation  of  at  a low- 
coordination  transition-metal  site  apply  equally  well  for  the  dissociation 
of  D2  at  such  a site. 

The  above  described  electronic  structure  of  the  L2MH2  (or  L^MD^) 
coordination  complex  leads  to  possible  explanations  of  the  observed  homo- 
geneous and  heterogeneous  catalytic  reactivity  of  (or  D^).  For  example, 
the  near  cancellation  of  the  contributions  of  the  bonding  (la^,  2a^)  orbitals 
and  antibonding  (3a^,  4a-j ) orbitals  to  metal -hydrogen  bond  strength,  leaving 
the  dissociative  2b2  bonding  orbital  dominant,  explains  the  relatively  weak, 
reversible  binding  of  H2  (or  D2)  to  such  complexes  and  their  ability  to 
activate  H2-D2  exchange  [13].  Since  such  a complex  is  also  a good  model  for 
H2  (or  D2)  dissociation  at  the  corner  atoms  of  a platinum  surface  step,  the 
results  suggest  why  atomic  steps  on  platinum  surfaces  are  essential  in 
dissociating  H2  and  D2  and  in  activating  H2~D2  exchange  [14]. 

The  electronic  structure  of  the  L2MH2  complex  also  suggests  a possible 
reaction  path  for  the  hydrogenation  of  unsaturated  hydrocarbons  at  low- 
coordination  transition-metal  sites.  The  4a^  orbital,  which  defines  the 
Fermi  energy  of  the  site,  is  closely  matched  in  symmetry,  energy,  and  electro- 
negativity to  the  tt  orbitals  of  hydrocarbons  such  as  acetylene  (C^)  and 
ethylene  (C2H4).  When  the  4a^  orbital,  which  is  an  antibonding  mixture  of 
L-M(dz2)  and  H2*  og  orbital  character,  is  only  partially  occupied  (as  is 
the  case  for  M = Ir,  Rh),  it  offers  a pathway  for  electron  flow  from  a C2H2 
(or  C2H4)  tt  orbital  to  the  dissociatively  "chemisorbed"  hydrogen.  Electron 
flow  directly  between  C2H2  (or  C2H4)  and  H2  in  the  gas  phase  via  the  filled 
tt  and  Og  orbitals  is  forbidden  by  the  Pauli  exclusion  principle,  whereas 
electron  flow  directly  between  the  tt  orbital  and  empty  au  orbital  is 
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forbidden  by  orbital  symmetry  [11].  Because  the  4a1  orbital  of  L2MH2  is 
antibonding  between  the  L^M  site  and  H2,  while  bonding  between  H2  and  C2H2 
(or  C2H^),  the  net  result  of  electron  flow  between  a n orbital  and  the  4a^ 
orbital  is  the  breaking  of  a C-C  tt  bond,  the  formation  of  two  new  C-H  bonds, 
and  the  expulsion  of  the  hydrogenated  species  C2H4  (or  C2Hg),  as  suggested 
by  the  reaction  path  shown  in  Fig.  3.  Also  shown  is  a contour  map  for  the 
4a1  orbital  of  the  L2MH2C2H2  reaction  intermediate  (the  third  step  of  the 
proposed  reaction  path)  formed  as  a result  of  the  interaction  of  acetylene 
with  the  L2MH2  complex.  The  incipient  formation  of  C-H  bonds  via  the  overlap 
of  the  C-C  tt  orbital  with  the  antibonding  metal -di hydride  orbital  and  the 
resulting  ethylene-like  configuration  are  clearly  visible  in  this  map.  It 
is  important  to  note  that  the  concerted  reaction  path  indicated  in  Fig.  3 
is  not  the  conventional  one  for  hydrogenation  on  ideal  transition-metal 
surfaces,  where  it  is  usually  assumed  that  chemisorption  of  acetylene  or 
ethylene  on  one  or  two  metal  sites  is  the  precursor  to  combining  with 
hydrogen  chemisorbed  on  neighboring  sites.  Nonconcerted  reaction  paths  in 
which  both  reactants  are  coordinated  to  the  same  metal  site  are  also  possible 
and  indeed  have  been  argued  to  be  favored  kinetically  in  certain  homogeneous 
reactions  [18].  Alternative  reaction  paths  at  low-coordination  transition- 
metal  sites  are  currently  under  investigation  in  conjunction  with  theoretical 
studies  of  the  reactivity  of  IrCl (CO) (Ph^P^  (Vaska's  complex)  [19]. 

In  this  section,  we  have  attempted  to  show  that  a detailed  theoretical 
study  of  the  electronic  structure  of  well  characterized  coordinatively 
unsaturated  transition-metal  complexes  and  their  interactions  with  H2  can 
not  only  lead  to  an  understanding  of  their  homogeneous  reactivity  but  can 
also  serve  as  a model  for  the  dissociative  chemisorption  and  heterogeneous 
reactivity  of  H2  on  low-coordination  transition-metal  surface  sites,  where 
definitive  structural  information  is  lacking.  There  are  many  useful  analogies 
to  be  made  between  molecular  transition-metal  coordination  complexes  and 
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surface-adsorbate  interactions  [20,21].  Such  analogies  are  probably  not 
fortuitous.  They  should  be  sought  after  and  the  common  basis  of  under- 
standing elucidated. 

In  concluding  this  section,  it  is  important  to  compare  the  present 
theoretical  approach  to  surface  reactivity,  based  on  molecular-orbital 
indices,  with  other  theoretical  approaches  to  this  problem  and  to  chemical 
reactivity  in  general.  In  the  applications  of  traditional  methods  of  quantum 
chemistry  (e.g.,  Hartree-Fock,  Configuration-Interaction,  Valence-Bond,  etc.) 
to  reaction  kinetics  and  thermochemistry,  one  usually  focuses  directly  on 
the  calculation  of  total  energies,  total-energy  differences,  and  "potential 
surfaces"  for  the  reactants.  This  is  generally  a computationally  difficult 
and  costly  process  to  carry  out  over  the  various  possible  reaction  paths, 
even  for  the  simplest  reactions.  To  appreciate  the  magnitude  of  this 
problem,  one  need  only  consider  the  recent  status  of  the  quantitative  first- 
principles  determination  of  the  kinetics  of  one  of  the  simplest  gas-phase 
chemical  reactions,  namely,  hydrogen-deuterium  exchange,  ^ + D HD  + H 
[22].  This  situation  hardly  makes  one  confident  in  the  efficacy  of  quantum 
theory  to  "predict"  the  path  or  kinetics  of  surface-activated  reactions  on 
the  basis  of  potential -surface  computations.  Even  if  one  carries  out  total 
energy  calculations  for  only  a few  representative  molecular  configurations, 
rather  than  for  the  entire  potential  surface,  there  is  still  the  uncertainty 
associated  with  the  direct  subtraction  of  two  total  energies,  which  are 
usually  large  numbers,  to  obtain  a relatively  small  energy  difference  of 
chemical  significance. 

The  theoretical  approach  described  in  this  paper,  while  not  a substi- 
tute for  ajD  initio  potential-surface  calculations,  circumvents  many  of  the 
difficulties  associated  with  the  latter  approach.  By  placing  emphasis  on 
the  determination  of  molecular-orbital  indices  of  reactivity,  rather  than 
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total  energies,  for  realistic  transi tion-metal  coordination  complexes  and 
clusters  simulating  local  bulk  and  surface  configurations,  one  retains  the 
molecular-orbital  picture  which  chemists  have  traditionally  used  and  makes 
contact  with  the  band-structure  concept  of  solid-state  and  surface  physics. 

IV.  THE  INTERACTION  OF  ATOMIC  HYDROGEN  WITH  TRANSITION-METAL  INTERFACES 

The  interaction  of  hydrogen  with  Group-VIII  transition  metals  such  as 
Ni , Pd,  and  Pt  is  of  fundamental  importance  in  the  understanding  of  (1)  the 
dissociative  chemisorption  and  reactivity  of  hydrogen  on  the  surfaces  of 
these  metals  [23],  (2)  the  catalytic  activity  of  small  particles  and  clusters 
of  these  metals  [15],  and  (3)  the  solubility  of  atomic  hydrogen  in  these 
metals  [24].  In  a previous  paper  [4],  it  has  been  shown  that  the  electronic 
structures  of  small  globular  Cu,  Ni,  Pd,  and  Pt  clusters,  as  calculated  by 
the  SCF-Xa-SW  method,  exhibit  most  of  the  characteristics  of  the  corres- 
ponding bulk  crystalline  band  structures,  while  having  additional  features 
corresponding  to  the  finite  cluster  size  and  "surface"  atoms.  In  the  pre- 
ceding section,  it  has  been  shown  that  low-coordination  transition-metal 
sites,  such  as  those  of  a cluster  or  "stepped"  surface,  have  local  features 
of  electronic  structure  which  explain  why  such  sites  activate  H2  dissocia- 
tion and  catalyze  H2-D2  exchange  [14]. 

In  the  present  section,  we  describe  the  results  of  SCF-Xa-SW  studies 
for  the  electronic  structure  of  four-atom  tetrahedral  Ni,  Pd,  and  Pt  clusters 
containing  atomic  hydrogen,  carried  out  in  collaboration  with  R.  P.  Messmer, 
D.  R.  Salahub,  and  C.  Y.  Yang.  This  work  includes  the  first  application  of 
the  relativistic  Xa-SW  formalism  developed  by  Yang  and  Rabii  [25]  to  metal 
clusters.  The  cluster  configurations  chosen  for  study  have  the  advantage 
that  they  are  large  enough  to  represent  the  local  effects  on  electronic  band 
structure  of  embedding  dissociated  hydrogen  atoms  in  a surface  or  bulk 
interstitial  environment,  yet  small  enough  to  permit  the  resolution  of 
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individual  metal -hydrogen  bonding  orbitals.  As  will  be  shown,  the  results 
are  in  good  agreement  with  photoemission  spectra  for  hydrogen  chemisorbed 
on  the  (111)  surfaces  of  crystalline  Ni , Pd,  and  Pt  [17],  strongly  suggesting 
that  the  hydrogen  atoms,  as  a result  of  their  small  size,  might 
penetrate  the  (111)  surface  planes  and  become  embedded  in  the  tetrahedral 
interstices  of  the  substrate.  Moreover,  the  computed  electronic  structures, 
in  conjunction  with  the  concept  of  orbital  electronegativity,  suggest  that 
in  those  metals  where  hydrogen  solubility  is  high  [24],  absorption  and 
chemisorption  are  closely  linked  and  can  be  described  by  essentially  identical 
theoretical  models.  This  is  supported  by  the  observation  that  the  photo- 
emission spectrum  for  hydrogen  dissolved  in  bulk  palladium  [26]  is  very 
similar  to  that  for  hydrogen  chemisorbed  on  palladium  [17]. 

Molecular-orbital  calculations  were  carried  out  for  the  representative 
nickel,  palladium,  and  platinum  clusters  using  both  the  standard  nonrela- 
tivistic  version  of  the  SCF-Xa-SW  method  and  the  recently  developed  rela- 
tivistic version  [25],  constraining  metal-metal  internuclear  distances  to 
be  equal  to  those  for  the  corresponding  crystalline  metals.  The  resulting 
relativistic  orbital  energies  for  the  tetrahedral  clusters  with  and  without 
interstitial  atomic  hydrogen  are  shown  in  Fig.  4.  Also  shown,  for  comparison, 
is  the  SCF-Xa  ls-orbital  energy  for  the  isolated  hydrogen  atom.  The  elec- 
tronic energy  levels  are  labeled  according  to  the  principal  partial-wave 
(s,  p,  d)  character  of  the  associated  molecular  orbitals,  and  the  highest 
occupied  orbital  in  each  cluster  is  indicated  by  the  "Fermi  level"  (ep). 

Since  these  clusters  are  intended  to  simulate  the  local  interstitial  bonding 
configurations  of  isolated  hydrogen  atoms  embedded  in  an  otherwise  perfect 
bulk  or  surface  lattice,  the  energy  levels  of  the  clusters  containing  hydrogen 
have  been  shifted  with  respect  to  those  of  the  corresponding  hydrogen- 
free  clusters  so  that  the  respective  Fermi  levels  line  up.  This  approxima- 
tion is  based  on  the  assumption  that  hydrogen  chemisorption  or  absorption 
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in  the  dilute  limit  does  not  severely  perturb  the  chemical  potential  of 
the  metallic  host  and  is  supported  by  the  observation  from  photoemission 
data  that  the  work  functions  of  crystalline  nickel,  palladium,  and  platinum 
change  by  no  more  than  + 0.2  eV  with  hydrogen  chemisorption  [17]. 

The  electronic  structure  of  each  metal  cluster  shown  in  Fig.  4 is 
characterized  by  a manifold  of  closely  spaced  d-levels  bracketed  by  s,p- 
or  s,p,d-hybrid  levels.  This  is  similar  to  the  results  obtained  for  larger 
Ni , Pd,  and  Pt  clusters,  as  described  in  Ref.  [4],  and  analogous  to  the 
overlap  of  the  "d-band"  by  the  "s,p-like  conduction  band"  in  the  re- 
spective bulk  crystal  line  metals  [27].  In  each  cluster,  the  Fermi  level  (ep) 
passes  through  the  top  of  the  d-band,  just  as  in  the  bulk  transition  metals 
[27].  Although  the  calculated  d-band  width  of  each  metal  cluster  is  less 
than  that  of  the  corresponding  bulk  metal,  the  trend  of  increasing  band- 
width and  downward  trend  of  the  energy  levels  from  Ni^  to  Pd4  to  Pt^  is 
similar  to  the  trends  for  the  crystalline  metals.  Furthermore, 
the  electronic  structures  of  the  Pd^  and  Pt^  clusters  are  more  nearly  alike 
than  those  of  the  Ni^  and  Pd4  or  Ni^  and  Pt^  clusters,  consistent  with  the 
band  structures  of  the  corresponding  crystalline  metals.  The  deepest  energy 
levels  shown  in  Fig.  4 for  Pd^  and  Pt^  respectively,  associated  with  cluster 
orbitals  having  the  a^  representation  of  the  Trf  point  group  and  corresponding 
roughly  to  the  Bloch  band-structure  states  having  the  r1  representation  of 
the  crystal  space  group,  are  predominantly  d-like  with  a small  amount  of  s,p- 
hybridization.  In  contrast,  the  deepest  energy  level  for  the  Ni^  cluster 
shown  in  Fig.  4,  also  associated  with  an  a^  molecular  orbital  and  Bloch 
state,  is  predominantly  s-like,  but  with  significant  d-orbital  hybridization 
and  some  p-like  character.  These  differences  between  the  electronic  structure 
of  nickel  aggregates  and  the  electronic  structures  of  palladium  and  platinum 
aggregates  are  crucial  to  understanding  the  differences  in  the  photoemission 
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spectra  for  hydrogen  chemisorbed  on  these  metals  [17],  as  well  as  the  differ- 
ences among  these  metals  with  respect  to  hydrogen  solubility  [24]  and 
catalytic  reactivity  [23]. 

According  to  the  discussion  presented  in  Section  II  of  this  paper,  the 
relative  positions  of  the  SCF-Xa  orbital  energies  for  the  Ni^,  Pd^,  and  Pt^ 
clusters  with  respect  to  the  SCF-Xa  Is-orbital  energy  for  atomic  hydrogen, 
as  shown  in  Fig.  4,  are  a measure  of  the  relative  orbital-electronegativity 
and  chemical -potential  differences  between  these  metal  aggregates  and 
hydrogen.  Thus  the  covalent  bonding  of  atomic  hydrogen  at  the  cluster 
interstices  is  governed  principally  by  the  proximity  in  energy  (or  electro- 
negativity) and  concomitant  overlap  of  the  symmetry-conserving  [11]  a^ 
orbitals  near  the  bottom  of  the  Ni^,  Pd^,  and  Pt^  d-bands  with  the  H Is 
orbital.  The  main  result  is  the  splitting  off  of  a hydrogen -metal  bonding 
energy  level  of  a^  orbital  symmetry  from  the  bottom  of  the  d-band  of  each 
cluster,  accompanied  by  much  smaller  level  shifts  within  the  d-manifolds, 
as  indicated  in  Fig.  4 by  the  orbital  energies  for  the  Ni^H,  Pd^H,  and  Pt^H 
clusters  and  the  connecting  dashed  lines.  The  metallic  4s-like  component 
of  this  a1  orbital  is  largely  responsible  for  the  bonding  of  hydrogen  to 
the  nickel  aggregate,  as  indicated  by  the  partial-wave  decomposition  of 
the  orbital.  However,  the  contribution  of  the  3d-like  component  to  the 
bonding  is  not  negligible,  amounting  to  35%  of  the  Ni^-H  a^  bonding  orbital 
charge.  This  result  is  inconsistent  with  the  claims  of  other  workers  [28,29] 
who,  on  the  basis  of  theoretical  studies  of  the  interaction  of  hydrogen  with 
only  one  or  two  nickel  atoms,  find  that  the  Ni  3d  orbitals  remain  essentially 
localized  and  atomic-like  and  therefore  conclude  that  these  orbitals  do  not 
contribute  to  the  chemisorption  of  hydrogen  on  nickel.  In  contrast  to  the 
results  for  nickel,  the  metal  d-orbital  components  almost  exclusively  dominate 
the  bonding  of  hydrogen  to  palladium  and  platinum  aggregates,  the  contributions 
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of  Pd  5s-  and  Pt  6s -1 i ke  components  amounting  to  only  15%  and  12%  for  Pd^H 
and  Pt^H  respectively.  The  covalent  overlap  of  directed  Pd  d-orbitals  with 
the  H Is  orbital  is  implicit  in  the  contour  map  of  the  a1  bonding  orbital 
shown  in  Fig.  5,  plotted  in  the  plane  defined  by  the  hydrogen  atom  and  two 
palladium  atoms.  These  results  underscore  the  danger  of  making  general 
conclusions  about  the  dominance  of  s-orbitals  over  d-orbitals  in  determining 
the  chemisorption  and  catalytic  reactivity  of  hydrogen  on  transition  metals 
exclusively  on  the  basis  of  theoretical  studies  of  first-row  transition 
metals,  as  has  recently  been  done  by  some  workers  in  the  published  litera- 
ture [28,29].  The  present  findings  are  essentially  unaltered  for  hydrogen 
interacting  with  nickel,  palladium,  and  platinum  aggregates  having  other 
cluster  configurations  large  enough  to  realistically  simulate  the  bulk  and 
surface  electronic  structures  of  these  metals  and  small  particles  thereof. 

For  example,  the  partial-wave  decompositions  of  the  orbitals  responsible 
for  the  bonding  of  hydrogen  at  the  octahedral  interstices  of  six-atom  clusters 
are  very  similar  to  the  results  described  above. 

The  most  striking  confirmation  of  these  theoretical  results  is  the 
photoemission  spectra  recently  measured  by  Demuth  [17]  for  hydrogen  chemi- 
sorbed on  the  (111)  faces  of  nickel,  palladium,  and  platinum.  For  each  metal, 
the  data  clearly  show  a chemisorption-induced  photoemission  peak  at  an  energy 
slightly  higher  than  the  metal  d-band  photoemission  peaks,  suggestive  of  a 
hydrogen-metal  bonding  state  or  "resonance"  split  off  from  the  manifold  of 
d-orbitals  as  predicted  in  Fig.  4.  On  the  basis  of  the  intensity  and  width  of 
the  chemisorption-induced  photoemission  peak  as  a function  of  incident  photon 
energy,  Demuth  [17]  concludes  that  hydrogen  chemisorption  on  Ni(lll)  occurs 
primarily  via  the  s-orbitals  (with  some  d-orbital  participation),  whereas  the 
metal  d-orbitals  dominate  hydrogen  chemisorption  on  Pd (111)  and  Pt(lll).  This 
interpretation  is  completely  consistent  with  the  partial-wave  decomposition  of 
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the  hydrogen-metal  bonding  orbitals  of  a1  symmetry  described  above.  Further- 
more, the  chemisorption-induced  spectral  modification  of  d-band  photoemission 
from  Pd(lll)  and  Pt(lll)  surfaces  [17]  is  explicable  in  terms  of  the  calculated 
shifts  of  d-manifold  energy  levels  in  going  from  Pd^  to  Pd^H  and  from  Pt^  to 
Pt^H  (Fig.  4).  On  the  other  hand,  the  uniform  enhancement  of  the  d-band  photo- 
emission from  the  Ni(lll)  surface  upon  hydrogen  chemisorption  is  consistent 
with  the  almost  negligible  shifts  calculated  for  the  d-orbital  manifold  in 
going  from  Ni^  to  Ni^H.  The  relative  energies  of  the  metal  cluster  d-orbital 
manifolds  with  respect  to  the  H Is  level  indicated  in  Fig.  4,  in  conjunction 
with  the  concept  of  orbital  electronegativity,  imply  that  a platinum  aggregate 
is  somewhat  of  an  electron  acceptor  with  respect  to  atomic  hydrogen,  whereas 
nickel  and  palladium  aggregates  are  slight  electron  donors.  This  finding  is 
consistent  with  Demuth's  [17]  photoemission  results  for  hydrogen  chemisorption 
on  the  surfaces  of  these  metals,  which  show  a slight  decrease  in  the  work 
function  of  Pt(lll)  and  slight  increases  in  the  work  functions  of  Ni(lll)  and 
Pd(lll)  upon  hydrogen  chemisorption.  The  close  correspondence  between  the 
theoretical  results  for  hydrogen  bonded  to  the  interstices  of  tetrahedral  Ni^, 
Pd^,  and  Pt^  clusters  and  the  photoemission  spectra  for  hydrogen  chemisorbed 
on  Ni(lll),  Pd(lll),  and  Pt(lll)  strongly  supports  our  conjecture  that  such 
chemisorption  leads  to  incorporation  of  hydrogen  atoms  in  the  tetrahedral 
interstices  bounded  by  the  (111)  surfaces. 

In  regard  to  the  well  known  fact  that  atomic  hydrogen  is  more  soluble 
in  palladium  than  in  nickel  or  platinum  [24],  the  almost  perfect  tuning  of 
the  palladium  cluster  d-orbital  electronegativities  to  the  hydrogen  ls-orbital 
electronegativity,  as  indicated  in  Fig.  4 by  the  relative  positions  of  the 
corresponding  energy  levels,  suggests  almost  perfect  covalency  between 
palladium,  in  aggregate  form,  and  atomic  hydrogen.  In  contrast,  nickel  and 
platinum  aggregates  are  respectively  electropositive  and  electronegative 
with  respect  to  hydrogen.  The  strength  of  a heteronuclear  chemical  bond. 


-21- 


as  originally  described  by  Pauling  [30],  can  be  viewed  as  having  covalent 
and  ionic  contributions  in  general.  It  is  well  known  that  the  solubility 
of  an  impurity  in  a metal  generally  decreases  with  increasing  electronegativity 
difference  between  solute  and  solvent,  other  factors  such  as  atomic 
size  remaining  constant  [31].  Thus  the  attainment  of  nearly  zero  net  orbital 
electronegativity  difference  between  palladium  aggregates  and  atomic  hydrogen, 
thereby  minimizing  ionic  contributions  to  the  bonding  and  optimizing  Pd( 4d) - 
H(ls)  covalency,  is  consistent  with  the  higher  solubility  of  hydrogen  in 
palladium,  as  compared  with  nickel  and  platinum.  The  labile  exchange  of 
dissociatively  chemisorbed  hydrogen  atoms  between  the  surface  and  underlying 
interstices,  making  the  metal  a reservoir  for  atomic  hydrogen,  could  facilitate 
the  reactivity  of  hydrogen  with  other  chemisorbed  molecules,  offering  a possible 
explanation  of  why  palladium  is  an  order  of  magnitude  more  active  in  catalyzing 
hydrogenation  reactions  than  nickel  or  platinum  [32].  Although  the  metal-metal 
internuclear  distances  in  the  clusters  have  been  constrained  in  the  present 
studies  to  the  values  for  the  corresponding  bulk  crystalline  metals,  previous 
theoretical  work  on  the  cohesive  energies  of  metal  clusters  has  shown  that 
the  equilibrium  internuclear  distance  decreases  somewhat  with  decreasing 
number  of  atoms  in  the  cluster  [33].  The  latter  finding  suggests  that  the  size 
factor  for  small  metal  clusters  is  less  favorable  for  Interstitial  hydrogen 
incorporation  than  larger  particles  or  crystallites,  thus  providing  a possible 
explanation  of  the  observed  reduction  of  hydrogen  solubility  with  decreasing 
particle  size  [34]. 

Finally,  the  SCF-Xa-SW  values  for  the  relative  4s-  and  3d-like  partial  - 
wave  components  of  the  Ni^-H  a^  bonding  orbital  have  already  been  used  by 
Schonhammer  [35]  to  parameterize  an  Anderson-type  Hamiltonian  for  hydrogen 
chemisorption  on  nickel.  From  this,  Schonhammer  has  calculated  the  adsorbate 
Green's  function  and  the  width  of  the  chemisorption-induced  photoemission 
peak,  yielding  a result  in  excellent  agreement  with  the  measurements  of  Demuth  [17]. 
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V.  THE  SURFACE  REACTIVITY  OF  IRON 

In  the  preceding  two  sections,  we  have  attempted  to  show,  applying  the 
concept  of  orbital  electronegativity,  how  local  coordination  and  composition 
influence  reactivity  at  transition-metal  interfaces.  In  the  present  section, 
by  extending  orbital  electronegativity  to  spin  orbitals  as  described  in 
Section  II,  we  wish  to  show  that  local  magnetic  spin  polarization  can  also 
affect  reactivity  at  transition-metal  interfaces,  using  the  surface  reactivity 
of  representative  iron  clusters  as  an  illustrative  example.  Of  all  the 
transition  metals,  iron  is  one  of  the  most  reactive.  In  its  pure  crystalline 
form,  iron  is  readily  oxidized.  Among  the  first-row  transition  metals,  iron 
is  the  most  active  one  for  dissociating  molecular  nitrogen  and  is  widely  used 
as  a catalyst  for  the  synthesis  of  ammonia  [36,37].  Iron  is  also  an  excellent 
catalyst  for  the  Fischer-Tropsch  synthesis  of  hydrocarbons  from  carbon 
monoxide  [23]. 

Spin-unrestricted  SCF-Xa  calculations  have  been  carried  out  by  C.  Y. 

Yang  [38]  for  the  electronic  structures  of  4-,  8-,  9-,  13-,  and  15-atom 
clusters  of  iron  having  tetrahedral  and  cubic  geometries.  All  clusters  show 
significant  spin  polarization  of  the  orbitals,  with  the  9-  and  15-atom 
clusters  of  bcc  geometry  exhibiting  spin-polarized  electronic  structures  that 
are  remarkably  similar,  to  the  extent  that  such  comparisons  can  be  made,  to 
the  ferromagnetic  band  structure  of  bulk  crystalline  a-iron  calculated  by 
Tawil  and  Callaway  [39].  For  the  basis  and  justification  of  comparisons 
between  cluster  electronic  structures  and  crystalline  band  structures,  which 
have  previously  been  made  for  copper,  nickel,  palladium,  and  platinum,  the 
reader  should  consult  Ref.  [4].  The  SCF-Xa  orbital  energies  for  Fe<j  and 
Fej5,  labeled  according  to  the  irreducible  representations  of  the  0h  symmetry 
group,  are  shown  in  Figs.  6 and  7,  respectively.  The  Fermi  level  is  indicated 
by  an  arrow.  The  differences  between  the  results  for  Fe-^  and  Fe^  are  mainly 
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quantitative,  rather  than  qualitative,  both  clusters  having  the  following 
properties  which  are  directly  comparable  with  the  band  structures  and  measured 
physical  properties  of  ferromagnetic  a-iron: 

(1)  The  Fermi  level  of  each  cluster  passes  through  the  center  of  the  manifold 
of  minority-spin  d-levels,  which  is  similar  to  the  intersection  of  the 
Fermi  level  with  the  minority-spin  d-band  of  ferromagnetic  a-iron  [39]. 

(2)  The  combined  density  of  levels  for  each  cluster  is  qualitatively  similar 
to  the  density  of  states  for  a-iron  deduced  from  photoemission  spectra 

[40] .  In  particular,  the  density  of  states  for  Fe^5  is  in  better  quanti- 
tative agreement  with  the  photoemission  data  [40]  in  regard  to  bandwidth 
than  is  the  density  of  states  derived  from  band-structure  calculations  [39] 
for  crystalline  iron,  suggesting  that  the  cluster  spin  orbitals  are  a 

more  realistic  description  of  the  initial  states  of  photoemitted  electrons 
than  are  delocalized  Bloch  states. 

(3)  The  magneton  numbers  per  atom  are  2.9  and  2.5  for  Feq  and  Fe1(-. 
respectively,  suggesting  convergence  to  the  2.2  value  characteristic 
of  ferromagnetic  a-iron.  The  SCF-Xa  results  for  successively  smaller 
iron  clusters,  such  as  Feg  and  Fe^,  indicate  convergence  to  larger  values 
of  the  magneton  number  approaching  the  atomic  limit. 

(4)  The  partial -wave  decomposition  of  the  Feg  and  Fe^g  cluster  spin  orbitals 
indicates  that  the  contribution  of  the  4s-like  orbitals  to  spin  polariza- 
tion, although  relatively  small,  is  opposite  in  direction  to  the  3d-like 
contribution,  in  good  agreement  with  neutron  diffraction  measurements 

[41] . 

(5)  The  spin  densities  along  the  [100]  directions  of  the  Feg  and  Fe^g  clusters, 
corresponding  to  the  direction  of  easy  magnetization  in  crystalline  a- 
iron,  are  significantly  larger  than  the  spin  densities  along  the  [111] 
directions,  due  to  the  greater  concentration  of  spin  density  in  the  eg 
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orbitals  than  in  the  t2g  orbitals.  This  result  is  in  excellent  quantita- 
tive agreement  with  neutron  diffraction  data  [41]  and  provides  a "real- 
space”  interpretation  of  the  latter. 

(6)  The  transition  from  ferromagnetism  to  paramagnetism  in  crystalline  a- 
iron  around  the  Curie  temperature  can  be  explained  within  the  framework 
of  the  Fe-j 5 cluster  model  in  terms  of  thermally  induced  long-range 
disordering  of  spin  clusters  and  localized  excitations  of  electrons 
within  each  cluster  from  the  "spin-down"  orbitals  at  the  Fermi  energy 
to  the  unoccupied  7t2gt  and  6egf  orbitals  lying  just  above  the  Fermi 
energy  (see  Fig.  7).  The  small  increase  of  net  magnetic  moment  of  each 
cluster  arising  from  the  thermally  induced  depletion  of  minority-spin 
orbitals  and  population  of  majority-spin  orbitals  is  consistent  with 
neutron  diffraction  measurements  in  the  vicinity  of  the  Curie  temperature. 
This  is  the  first  such  explanation  of  a magnetic  phase  transition  and 
lies  beyond  the  scopes  of  conventional  band  theory  and  molecular-field 
theory. 

A more  detailed  discussion  of  the  relationships  between  these  cluster 
calculations  and  the  physical  properties  of  crystalline  a-iron  will  be  the 
subject  of  a forthcoming  paper  [42],  including  a study  of  the  effects  of  differ- 
ent boundary  conditions  on  the  clusters,  a comparison  of  cluster  and  bulk 
densities  of  states,  and  calculations  of  the  Heisenberg  exchange  integral  and 
Curie  temperature  via  Slater's  [6]  transition-state  description  of  localized 
spin  excitations.  It  is  important  here  to  underscore  the  fact  that,  despite 
the  finite  molecular  nature  of  the  clusters  and  the  appearance  of  certain 
"surface-related"  features  of  the  cluster  electronic  structures  (see  below), 
the  spin-polarized  SCF-Xa  results  for  the  bcc  Fe^  and  Fe^  clusters  provide 
a remarkably  successful  model  for  the  electronic  and  magnetic  properties  of 
crystalline  a-iron.  Similar  theoretical  studies  have  been  made  for  cubo- 
octahedral  Fe13  clusters  representing  fee  crystalline  y-iron  [38].  The  above 
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results,  along  with  SCF-Xa  cluster  models  (containing  up  to  44  atoms)  for 
crystalline  aluminum  and  chemisorption  thereon  [43],  contradict  the  work  of 
van  Dyke  [44]  who  concludes  that  such  cluster  calculations  do  not  yield  an 
adequately  converged  description  of  the  properties  of  crystalline  metals. 

We  can  benefit  by  the  large  effective  surfaces  presented  by  the  Feg  and 
Fe-^  clusters  to  discuss  the  surface  reactivity  of  iron,  as  has  already  been 
done  for  nickel,  palladium,  and  platinum  in  Ref.  [4]  and  in  the  preceding 
section  of  this  paper.  Among  the  manifolds  of  densely  spaced  d-orbital  energies 
for  the  Feg  and  Fe^  clusters  (see  Figs.  6 and  7)  are  levels  which  correspond 
to  antibonding  spin  orbitals  primarily  localized  on  and  spatially  oriented 
away  from  the  cluster  boundaries  or  "surfaces,"  especially  in  the  vicinity  of 
the  Fermi  energy  in  the  minority-spin  d-manifold.  Many  of  these  spin  orbitals 
have  the  proper  spatial  character  for  symmetry-conserving  [11]  overlap  with 
the  orbitals  of  certain  reactant  molecules.  Moreover,  the  spin  polarization 
raises  the  minority-spin  orbitals  to  higher  energies  in  comparison  with  the 
non-spin-polarized  limit,  effectively  reducing  the  orbital  electronegativity 
and  facilitating  overlap  of  the  highest  occupied  iron  surface  spin  orbitals 
with  the  lowest  unoccupied  or  partially  occupied  orbitals  of  reactant  molecules 
such  as  N^,  CO,  and  O2.  Since  the  latter  orbitals  are  antibonding,  overlap 
and  effective  electron  flow  from  the  iron  surface  spin  orbitals  to  the 
unfilled  reactant  orbitals  should  promote  molecular  dissociation,  the  precursor 
to  surface  reactions  of  these  molecules,  such  as  ammonia  synthesis,  Fischer- 
Tropsch  synthesis,  and  surface  oxidation  [23].  This  argument  is  clarified  by 
the  direct  comparison  of  representative  Feg  cluster  spin-orbital  energies  with 
the  orbital  energies  of  Ng,  CO,  and  0^  in  Fig.  8,  which  is  equivalent  to 
comparing  orbital  electronegativities  as  discussed  in  Section  II.  Also  included 
for  comparison  are  the  SCF-Xa  orbital  energies  of  a Pt^  cluster  previously 
shown  to  exhibit  many  of  the  characteristics  of  the  bulk  and  surface  electronic 
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structures  of  crystalline  platinum  or  small  particles  thereof  [4].  In 
contrast  to  iron,  there  is  no  magnetic  spin  polarization  in  the  platinum 
cluster,  corresponding  to  the  nonmagnetic  character  of  bulk  platinum,  and 
the  high  density  of  d-orbitals  around  the  Fermi  level  is  poorly  matched  in 
energy  and  orbital  electronegativity  to  the  lowest  unoccupied  N2  and  CO 
orbitals.  This  result  is  consistent  with  the  experimental  fact  that  platinum 
is  significantly  less  active  than  iron  in  promoting  N2  or  CO  dissociation  and 
in  catalyzing  ammonia  or  Fischer-Tropsch  synthesis  [15,23].  Nickel  is  also 
magnetic,  but  the  spin  polarization  of  the  d-orbitals  is  not  large  enough  to 
yield  a density  of  states  around  the  Fermi  level  which  is  as  well  matched  in 
energy  and  spin-orbital  electronegativity  to  the  lowest  unoccupied  orbitals 
of  N2  or  CO  as  in  the  case  of  iron.  This  comparison  is  made  in  Fig.  9 and  is 
consistent  with  the  fact  that  nickel  is  a poor  ammonia  synthesis  catalyst  and 
is  less  active  than  iron  in  promoting  Fischer-Tropsch  synthesis  [15,23].  With 
respect  to  the  latter  class  of  reactions,  nickel  is  used  mainly  as  a methana- 
tion  catalyst  [15].  The  fact  that  iron  is  more  readily  oxidized  than  platinum 
is  explained  by  the  relative  differences  of  the  corresponding  orbital  energies 
and  electronegativities  with  respect  to  the  partially  occupied  antibonding 
Tig  orbital  of  02,  as  shown  in  Fig.  8. 

In  an  attempt  to  elucidate  further  the  mechanism  of  ammonia  synthesis  on 
iron  surfaces,  Yang  [38]  has  constructed  cluster  models  for  the  so-called  C7 
active  site,  present  on  the  (111)  plane  of  bcc  iron  [37],  and  for  a surface 
nitride  of  the  type  which  may  be  formed  after  the  rate-limiting  step  of  N2 
dissociation.  In  Fig.  10,  the  Xa  spin-polarized  orbital  energies  of  an 
Fe^Ng  cluster  representing  a face-centered  surface  nitride  are  compared  with 
the  orbital  energies  of  the  pure  Fe^  cluster.  The  distinguishing  feature 
of  the  nitride  electronic  structure  is  a relatively  narrow  band  of  N 2p-like 
levels  near  the  bottom  of  the  Fe  d-band,  as  indicated  in  Fig.  10,  a result 
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which  is  confirmed  by  recently  measured  photoemission  spectra  for  nitrogen 
chemisorbed  on  iron  [40].  A detailed  comparison  of  the  FegNg  and  Feg 
electronic  structures,  along  with  the  measured  change  in  work  function  of 
iron  upon  nitrogen  chemisorption  [40],  suggests  some  electronic  charge 
transfer  from  iron  to  nitrogen  in  the  formation  of  the  surface  nitride.  The 
effective  negative  charge  on  the  surface  nitrogen  atoms,  together  with  the 
close  matching  of  the  nitride  and  H Is  orbital  energies  and  electro- 
negativities, should  facilitate  protonation  of  the  nitrogen  atoms,  the  formation 
of  N-H  bonds,  and  ultimately  ammonia  synthesis. 

A similar  argument  applied  to  the  interaction  of  CO  with  iron  leads  one 
to  conclude  that  a surface  iron  carbide,  formed  from  the  dissociation  of  CO 
on  the  iron  surface,  could  facilitate  protonation  of  the  surface  carbon 
atoms,  the  formation  of  C-H  bonds,  and  the  Fischer-Tropsch  synthesis  of 
hydrocarbons.  There  is  already  evidence  that  iron  carbide  will  catalyze 
Fischer-Tropsch  synthesis  [23].  A theoretical  study  of  this  surface  reaction 
is  in  progress  at  M.I.T.,  as  are  studies  of  the  reactivity  of  N2  on  ruthenium 
and  osmium,  the  second-  and  third-row  transition  metals  which  are  good 
catalysts  for  ammonia  synthesis. 
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Figure  3. 
Figure  4. 

Figure  5. 
Figure  6. 

Figure  7. 
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Figure  9. 

Figure  10. 


SCF-Xa  orbital  energies  for  coordinati vely  unsaturated  transition- 
metal  complexes  representing  low-coordination  transition-metal 
sites  and  dissociative  hydrogen  chemi sorption  thereon.  The  highest 
occupied  orbital  is  indicated  by  the  "Fermi  level"  Cp. 

Contour  maps  of  the  principal  bonding  and  antibonding  molecular- 
orbital  wavefunctions  corresponding  to  the  orbital  energies  of 
the  L2MH2  complex  shown  in  Fig.  1. 

Possible  reaction  path  for  the  hydrogenation  of  acetylene  at  a 
coordinati vely  unsaturated  transition-metal  site. 

Relativistic  SCF-Xa  orbital  energies  of  tetrahedral  Group-VIII 
transition-metal  clusters  with  and  without  interstitial  atomic 
hydrogen. 

Contour  map  for  the  a^  bonding  orbital  of  the  Pd^H  cluster. 
Spin-polarized  SCF-Xa  orbital  energies  of  a 9-atom  bcc  iron 
cluster. 

Spin-polarized  SCF-Xa  orbital  energies  of  a 15-atom  bcc  iron 
cluster. 

Comparison  of  the  SCF-Xa  orbital  energies  of  the  Feg  and  Pt13 
clusters  with  those  of  the  N2,  CO,  and  02  molecules. 

Comparison  of  the  SCF-Xa  orbital  energies  of  the  Feg  and  Ni^3 
clusters  with  those  of  the  N2  and  CO  molecules. 

Comparison  of  the  SCF-Xa  orbital  energies  of  an  FegNg  cluster 
representing  a "surface  iron  nitride"  with  the  orbital  energies 
of  an  Feg  cluster  representing  pure  iron  and  with  the  orbital 
energies  of  atomic  hydrogen  and  nitrogen. 
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